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Abstract 

Shape analysis methods have in the past few years become very pop¬ 
ular, both for theoretical exploration as well as from an application point 
of view. Originally developed for planar curves, these methods have been 
expanded to higher dimensional curves, surfaces, activities, character mo¬ 
tions and many other objects. 

In this paper, we develop a framework for shape analysis of curves in 
Lie groups for problems of computer animations. In particular, we will 
use these methods to find cyclic approximations of non-cyclic character 
animations and interpolate between existing animations to generate new 
ones. 
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1 Introduction 

Motions of virtual characters in movies and interactive applications are usu¬ 
ally represented using a skeletal animation approach where the data consists of 
curves tracking the positions of the bones throughout the motion. These curves 
can be processed by mathematical methods to produce new motions [36]. It 
has previously been shown that shape analysis methods can be successfully 
applied to solve problems in computer animation by considering entire anima¬ 
tions as curves and shapes [6, 15]. In practice, the data consists of curves in 
SO(3) d , where d is the number of bones in the skeleton. However, in the ear¬ 
lier approaches, curves have been represented using Euler angles, neglecting the 
underlying Lie group structure. We here report on the results we obtained by 
appropriately including the underlying geometric structure in the mathematical 
models and their numerical discretizations. The intrinsic geometric formulation 
is robust and works very well in problems of motion blending and curve closing, 
where earlier the same performance could be only obtained by using ad hoc 
strategies, e.g. keeping track of carefully chosen feature points along the curves. 

In section 2, we will briefly introduce shape analysis and motivate how tech¬ 
niques from shape analysis on Lie groups can be applied to computer animation 
by treating character animations as points in an infinite dimensional manifold. 
This manifold is in fact an infinite dimensional Lie group where we are interested 
in computing distances and geodesics. 

In Section 3, we discuss some of the main tools for shape analysis on Lie 
groups, which are later applied in this paper. An approach to curves evolving 
on Riemannian manifolds was earlier presented in [44] . In our work, we exploit 
the additional structure provided by the Lie group setting. To perform anima¬ 
tion blending (i.e., interpolation between existing animations to create a new 
one), we are interested in simple and computationally efficient approaches to 
compute geodesics between two given shapes and to this end we define a metric 
on shape space. We show that this metric is associated to the geodesic distance 
of the pullback of an L 2 -inner product (cf. Theorem 3.11 and Theorem 3.16). 
The Lie group formalism allows us in Section 4 to efficiently solve the curve 
closing problem (which we use to approximate an existing non-periodic anima¬ 
tion with a periodic one) using a gradient flow approach. Finally, in Section 
5, we present numerical results both for problems in animation blending and 
animation closing. 
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We demonstrate that the proposed techniques exhibit better qualitative per¬ 
formance compared to previous work [6, 15]. In animation blending the new 
approach allows interpolating between a wider range of motions. In animation 
closing the proposed Lie group formulation is naturally intrinsic, and allows 
avoiding undesired artefacts due to coordinate singularities. While our specific 
applications use the special orthogonal group SO( 3) as the underlying Lie group, 
the techniques developed in this paper are not restricted to this setting. 

2 Background 

2.1 Shape analysis 

Many problems in object and activity recognition can be formulated in terms 
of similarities of shapes [3, 12, 15, 16, 23, 34, 39, 41, 42, 45, 46]. By shapes 
we typically mean unparametrized curves in a vector space or on a manifold, 
although similar methods have been developed for surfaces [2, 7, 29, 30]. 

In recent years, a number of methods based on differential geometry have 
been developed to tackle such problems. Of particular importance is the ques¬ 
tion of how to model and work with unparametrized curves. 

A popular approach is to define shapes as equivalence classes of certain map¬ 
pings, where the equivalence relation is induced by reparametrization. Given 
two curves c 0 ,ci : I —>• M with / = [a, b] C R and M a vector space or a 
manifold, we define equivalence classes [co], [ci] via the equivalence relation: 

c 0 ~ ci 3 ip : c 0 = Ci o ip, 

where <p is a smooth, strictly increasing bijection on I. We denote by V the 
space of parametrized curves containing co and Ci. Typical choices of V include 
absolutely continuous functions, immersions, embeddings or piecewise linear 
functions [4, 31]. The equivalence classes, or shapes, can then be collected in 
the corresponding shape space: 


S :=V/~. 

Shape analysis then concerns itself with the study of the spaces S and V 
from both a theoretical and a practical perspective. 

Many applications require for example a distance function in S to measure 
similarities between shapes, which can be used to perform statistical analysis 
such as clustering and object recognition [41]. Distance functions on shape space 
S are typically obtained from a distance function d-p on the underlying space V 
of parametrized curves as follows: 

ds(M,[ci]) := inf dp(co,Ci op), 

where p ranges over all possible curve reparametrizations. The computation 
of ds amounts to the solution of an optimization problem with appropriate 
numerical techniques [38, 42]. 

In recent years, a distance function on V based on the so-called elastic metric 
has become popular. The elastic metric is a first-order Sobolev-type metric, 
which is easy to compute and has desirable theoretical properties [4, 7, 33, 42], 
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In this paper, we consider curves and shapes in Lie groups to solve problems 
in computer animation. In that context the use of Lie groups arises naturally 
by using Euclidean transformations and rotations to describe motions of virtual 
characters. 

2.2 Skeletal animation 

In computer graphics, used and seen in movies, tv-series and video games, but 
also in educational and scientific software, virtual characters are most commonly 
represented as surfaces in M 3 . 

Motions of such characters are usually represented using a skeletal anima¬ 
tion approach. The underlying skeleton consists of bones connected by joints. 
The vertices of the surface mesh are attached to bones, i.e., their positions are 
specified in a coordinate system that is aligned with a bone. Then, when the 
skeleton is animated by specifying the positions of all the joints as a function 
of time, the vertices move accordingly. See Figure 1 for an example virtual 
character. 

A skeleton is a rooted tree, where each node (bone) represents a Euclidean 
coordinate system. The edge (joint) between a node and its parent represents 
a Euclidean transformation indicating how the two nodes’ coordinate systems 
are positioned and oriented relative to each other. By following this kinematic 
chain of relative transformations from a node to the root, all local coordinate 
systems can be brought into a single global coordinate system. 

A character’s pose is specified by assigning values to all degrees of freedom in 
the skeleton (i.e., all the joints). Each such configuration is an element of joint 
space J := SE(3) d . An animation is then a function from some time interval 
[a, b\ into joint space, specifying a character pose for every point in time. 

For human characters, bones have fixed lengths and joint space therefore consists 
only of rotations between bones (i.e., elements in 50(3), the special orthogonal 
group) instead of more general Euclidean transformations {SE{ 3)). 

A typical approach of generating such animations for use in videos or games 
is motion capturing , where an actor’s motions are recorded from multiple points 
of view and the underlying skeletal motion is extracted. This is often preferred 
over alternative means of generating animations, such as manual construction 
or inverse kinematics. 

Motion capturing however suffers from limitations inherent in its process, 
namely that the data is static and limited in range. If we want an animated 
character to run at different speeds, we need corresponding recordings. If we 
want a virtual character to keep walking forward for an indeterminate amount 
of time (e.g., in a video game under a player’s control), we need to find a way 
of adapting a finite walking animation so that we can repeat it without visible 
discontinuities. 

Much work has been done to procedurally manipulate motion data to tackle 
these and many more problems in computer graphics. We refer to, among many 
others, [19, 24, 25, 26, 36, 40] for an overview and some specific approaches in 
this field. 

To take into account the geometry inherent to skeletal animation, in the follow¬ 
ing sections, we propose methods for motion blending and animation closing on 
Lie groups. 
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Figure 1: Surface mesh depicting a human character with underlying skeleton. 
Only bones with a degree of freedom are shown, i.e., for example the hips are 
not visible in the skeleton and the legs appear at an offset from the character’s 
center. 
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3 Shape analysis on Lie Groups 

In this section we will develop a framework for shape analysis for curves on Lie 
groups. Albeit our main application is shape analysis for curves on SO( 3), the 
theoretical framework exhibited here is fairly general. The methods discussed 
in this and the next section can also be applied to certain classes of infinite¬ 
dimensional Lie groups. 

In the following, G will refer to a Lie group with Lie algebra (g, [•, •]), iden¬ 
tity e and multiplication • : G x G —> G. Here G might even be an infinite¬ 
dimensional Hilbert Lie group, i.e. a Lie group in the sense of [28] (or [35]) 
modelled on a Hilbert space. 1 We will denote left and right translations by L g 
and R g respectively, i.e., L gi (g 2 ) = 31-32 and R gi { 92 ) = 32 ■ 31 • We further 
denote by T the evolution operator , which is defined as 

T : C°°(/,g) ->{c£ C°°(I, G) : c(0) = e} =: C™ (I, G) 

T(q){t) ■■= c{t) where gjc(t) = R c (t)*(q(t)), c(0) = e, 

where I C R is an interval with 0 £ I and R gif = T e R g is the tangent map of the 
right translation at the identity e. Without loss of generality, we let I = [0,1]. 

The inverse of the evolution operator is the so called right logarithmic deriva¬ 
tive 


5 r : C™(I,G)^C°°(I, Q ), 

S r c := (R-'Uc). 

3.1 Manifolds of smooth mappings 

In this section we recall the construction of the manifold structure on spaces of 
smooth mappings with values in a Lie group. Moreover, we review some basic 
facts on these manifolds which will be used throughout the article. Following 
[27, Chapter 42], we construct the manifold structure on C°°(/, G ) using a local 
addition on G. A local addition allows us to choose local parametrisations on a 
manifold in a smooth way. 

Definition 3.1 (Local Addition, cf. [27, 42.4]). A local addition on a manifold 
M is a smooth map E: TM D fl —► M defined on an open neighbourhood fi of 
the zero-section such that 

• E(Cb) = x for all x £ M, where 0 X £ T X M is the zero element, 

• (7 r T M,S): n -» M x M, £ x >->• (7 r T M(£r)> s tez)) = (z, £(&:)) induces a 
diffeomorphism onto a neighbourhood of the diagonal in M x M. Here 
7 ttm is the tangent bundle projection. 

Since G is a Lie group we can construct a local addition as follows. Choose 
a chart ip: G D C g around the identity in G. Denoting elements in 

the tangent space over g £ G by £ g , we obtain a local addition via 

E: TG D n := |J R g ^W^) -> G, £ g 3 ■ ^“'(^fe))- 

geG 

1 Assuming that the Lie group is modelled on a Hilbert space will assure that our results 
carry over without any change. However, some of the results extend even to Banach Lie 
groups, cf. Remark 4.3. 
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Now define G£° (I, TG) := {/ £ G°°(/,TG) | tt tg o/ = a} for a £ G°°(/,G). 
This is a vector space with the pointwise operations (in the fibres T a ( t )G for 
each t). 

In the following we endow spaces of smooth functions G°°(J, M) to a (pos¬ 
sibly infinite-dimensional) manifold M with the compact open G°°-topology. 
This topology allows to control functions and their partial derivatives on any 
compact subset of /. We refer to [35, Section 1.5] for more information on this 
topology. Then U a := {h £ C°°(I,G) | Vf £ /, ( a(t),h(t )) £ Im(7TTGj51)} is 
open in C°°(I,G) and one can prove that the assignment 


V a : C°°(I, G)DU a -> C?(I,TG), / ( 7 T TG , S )- 1 o (a, /) 


is a manifold chart for C°°(I, G ). By [37, Theorem 7.8] (and the remarks before 
the cited theorem) these charts turn C°°(I,G) into a Frechet manifold (i.e. a 
manifold modelled on a locally convex space which is complete and metrizable). 

At this point we leave the realm of Banach manifolds, whence the standard 
definition for smooth maps (i.e. viewing the derivative as a continuous map 
to a space of continuous operators) breaks down. One way to define smooth 
maps beyond the Banach setting is the so called Bastiani calculus (see [1]): 
A map between Frechet spaces is smooth if all iterated directional derivatives 
exist and are continuous in a natural sense (see e.g. [35, 1.2] for more details). 
There are other (in general inequivalent) ways to define smooth mappings on 
infinite-dimensional spaces, such as the so called convenient calculus (see [27]). 
Fortunately, these choices yield the same smooth maps on Frechet manifolds. As 
we exclusively work with Frechet manifolds we can thus disregard the differences 
and freely use results formulated in both calculi . 2 

Furthermore, we quote from [27, 41.10, Theorem 42.13] and [37, Theorem 
7.9] the following results which we will use later on. 

Proposition 3.2. 

1. Pointwise multiplication and inversion induce a Frechet Lie group struc¬ 
ture on G°°(/,G). Its Lie algebra is C°°(I, 9 ) and we let ec- I —> G be 
its identity. 


2. The set Imm(/, G) is open in C°°(I,G ) and for an open subset U C G, 
the set C°°(I,U) := {/ £ G°°(I,G) \ f(I) C U} is open. 

3. By [17, 1.9], C%°(I,G) := {c £ G°°(J, G) : c(0) = e} is a closed Lie 
subgroup with Lie algebra {/ £ C°°(I,g) : /(0) = 0 }. 3 

4■ Using the bijection (-) A : G°°(K, G°°(/, G)) —> G°°(K x I, G), 7 a (e, t) := 
7 (e) (t) one identifies the tangent bundle TC°°(I,G) with C°°(I,TG). Fi- 
brewise this isomorphism is given by 



where 7 : K — >■ G°°(/, G) is smooth, with 7(0) = a. 

2 Note that we have already done this, as [37] (Bastiani calculus) just generalizes [27, 
Section 42] (convenient calculus). Further, both calculi can handle smooth maps on the non¬ 
open domain I = [a, 6], see e.g. [37, Definition 7.2] and [27, Chapter 24], 

3 Here we have used that ip ea = ■)/>*: C°°(I,V^) —> C°°(/, Mfo),/ 1 —> o /, where 
C^ChWfo) := {/ S C°°(/,fl) [ /(/) C W^,} is open in C°°(/, 0 ). The Lie group structure 
from 1. coincides with the one defined in [17]. 
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3.2 Shape space and distance functions 

We want to measure distances between two shapes, i.e., unparametrized curves 
on G. To do this, one represents unparametrized curves as equivalence classes 
of parametrized curves under certain reparametrizations. 

We will model parametrized curves in G as immersions (i.e., smooth func¬ 
tions with a nonvanishing first derivative) and denote the space of these curves 
by V := Imm(J, G). Then one can define the shape space S as the quotient 
space 

S:= P/Diff+(/), 

where Diff + (/) denotes the group of orientation preserving diffeomorphisms 
on /, acting on V from the right [4]. In the context of parametrized curves, 
Diff + (I) can be thought of as the group of all possible orientation preserving 
parametrizations of a curve. 

Our goal is to find a distance function ds on S. To this end, we consider an 
appropriate distance function on V . The distance of two parametrized curves 
will be measured as the infimum of the length of piecewise smooth curves con¬ 
necting these curves (i.e. the geodesic distance of a Riemannian metric). Note 
that in general, the distance will only be a pseudometric 4 , see Remark 3.5 be¬ 
low. However, the geodesic distance considered in our approach will turn out 
to be a metric (on suitable submanifolds). To assure that the distance descends 
to a distance function on S , the pseudometric needs to satisfies the following 
invariance property. 

Definition 3.3. Let dp : V x V —► [0, oo[ be a pseudometric. Then dp is called 
reparametrization invariant if 

dp(co,ci) = d-pico o <£>,ci o (p) \/ip g Diff + (J). (1) 

In other words dp is invariant with respect to the diagonal (right) action of 
Diff+(/) on V x V. 

Let [c 0 ], [ci] € S be equivalence classes and pick arbitrary representatives 
Co £ [co] and ci £ [ci]. This allows us to define a pseudometric ds as 

ds(M, [ci]) := inf d-p(co, C\ o <p). (2) 

V GDift+(7) 

Then we obtain the following result. 

Lemma 3.4. If dp is a reparametrization invariant pseudometric on V , then 
the pseudometric ds on S, as defined in Equation (2), is independent of the 
choice of representatives cq and c\. 

Proof. Due to the infimum over <p, the distance ds is independent of the choice 
of the representative C\. With respect to the choice of representative Co, we 

4 A pseudometric d: V X V —> R U {oo} satisfies all axioms of a metric but might fail to 
distinguish different points, i.e. in general d(x, y) = 0 will not imply x = y. 



easily find: 


<is([co ° i>\, [ci]) = inf d v (c 0 ot/>,ci o ip) 
veDi ff+(/) 

= inf d-pfco o ip o ip~ l , Ci o ip o tb~ ) 

(/) 

= inf dp(c 0 , ci 07 ) 

7£Diff+(/) 

= d 5 ([co], [ci]). □ 

Remark 3.5. Contrary to the finite-dimensional case, the geodesic distance of 
an infinite-dimensional Riemannian manifold might vanish. For example, it is 
well known that the geodesic distance with respect to the standard L 2 -metric on 
V vanishes everywhere. This is the reason why we formulated the above results 
using pseudometrics. 

However, it is also known that this pathology does not occur for the class of 
first order Sobolev metrics. We refer to [5] for more information and references. 

Next we will show how to obtain a reparametrization invariant geodesic 
distance function from a certain first order Sobolev metric. Thus its geodesic 
distance will turn out to be metric (on certain submanifolds of V). 


3.3 SRV transform for curves on a Lie group 

The main idea is to construct a well behaved mapping which allows us to pull 
the L 2 inner product on C°°(/, g) back to a Riemannian metric on a suitable 
submanifold of V. Inspired by the approach in [42, 43, 44], we define the map 


n : Imm(7, G) -> {q £ G°°(/,g) | q(t) ^ 0 Vt £ 1} = C°°(I,g\ {0}), 


q(t) = TZ{c){t) := 


V¥m 


( 3 ) 


where the norm || • || is induced by a right invariant metric on G. In particular, 
this entails that the norm on g is induced by an inner product (•, •) on g. Note 
that TZ is not injective since we lose information regarding the curve’s starting 
point. For G = this map is known as the Square Root Velocity Transform 
(SRVT) [42]. Hence we will also call TZ the square root velocity transform if V 
is a space of mappings with values in an arbitrary Hilbert Lie group. 

Let us first note some properties of the SRVT. 


Lemma 3.6. The SRVT 1Z is 

1. equivariant with respect to reparametrizations, i.e. V(c, if)ePx Diff + (I) 
we have 7 Z(c o (p) = lZ{c) o ip ■ i/ip, 

2. translation invariant, i.e. \/g £ G and c £ V we have 7Z(R g o c) = 7 Z(c). 
Proof. 1. A straight forward computation yields 


lZ(c o <p)(t) 


- R e (yw)«(c(y( f ))y , (*)) 




R cd(t )>(c(y>(*))) 

VWv<M\ 


2. This follows from R c( \ )g *{R g *c(t)) = R c ^(c{t)) (apply [27, 38.1 Lemma] 
to the constant map 1 1 —x g) and the definition of 7 Z. □ 
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Definition 3.7. Define a pseudometric on V = Imm(7, G) via 


dv(co, ci) 



qi{t)\\ 2 dt = d L 2 (TZ(c 0 ),TZ(ci)), 


( 4 ) 


where qi := 7Z(ci), i = 0,1. 

Notice that dp is only a pseudometric as it does not distinguish between 
cGV and R g o c (for g G G by Lemma 3.6). 

Proposition 3.8. The pseudometric dp from Definition 3.7 is reparametriza- 
tion invariant. 


Proof. In order to prove that dp is reparametrization invariant, we need to show 
that the property (1) holds. 

Using the definition (4) and the substitution s '■= <p(t), the reparametrization 
equivariance (see Lemma 3.6) implies: 


dp{c 0 oip, Cl o tp) = 


1 Viliam vwdi&m ) 


2 \ 1/2 
dt 


V- 


R c 0 (s)M 3 )) _ R C1 (S)M( S )) 


1/2 


ds 


= dp(c o, Ci). 


□ 


We would like to realize dp as the geodesic distance of a Riemannian metric 
on V which arises by pullback with TZ. In the classical case, of the abelian Lie 
group R 2 , the map TZ reduces to the SRVT considered in [3]. As was observed 
in [3, 3.4 and Remark 3.9], the map TZ is not even infinitesimally injective, i.e. 
the kernel of the tangent map T C 7Z is not trivial. Hence, pulling back the L 2 - 
metric will not result in a Riemannian metric on all of V. However, the map TZ 
restricts to a diffeomorphism onto a certain submanifold of V. Further, on this 
submanifold dp will coincide with the geodesic distance induced by the pullback 
of the L 2 -metric. Before we prove this, let us first introduce some auxiliary 
mappings and an explicit formula for the inverse of the SRVT. Consider the 
scaling maps 

sc: C™(/, g \ {0}) —> C™(I, g \ {0}), , ~ (. ~ (f) 

sc -1 : 0 \ { 0 }) ->• 0 \ { 0 }), (t ^ q(t)\\q(t)\\) 

Then it is easy to see that 7 Z(c) = sco<5 r (c) for all c G V. We will see in Lemma 
3.9 that TZ induces a diffeomorphism VC] Cf°(I,G) —> g\ {0}). By abuse 

of notation we write TZ~ 1 := (ft| PnC == ( / iG) ) _1 .Before we prove that TZ~ l is 
smooth and provide a formula, consider the set 

?» = {c £ Imm(7, G) : c(0) = e}. 
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Now V * = PnCf (/, G ) and by [17, 1.10] V * is a closed submanifold of V. Note 
that 7Z" 1 is the inverse of the SRVT when restricted to the submanifold V*. 

In the following, we restrict our investigation to 77*, i.e. we consider only 
curves starting at the identity element in G. We remark that this is only a mild 
restriction as the group operations in G allow us to transport any smooth curve 
c: I —> G to the smooth curve c(t) ■ c(0 ) _1 starting at the identity. 

Lemma 3.9. The scaling maps (5) are smooth diffeomorphisms. Moreover, 
the SRVT is smooth and induces a diffeomorphism V * — > C°°(I,q \ {0}) with 


TZ~ 1 : g \ {0}) C%° (/, G) 


dc 


K = r o sc 1 (q) = c(t), where —=R c{t> (q(t)\\q\\), c{V) = e. 


( 6 ) 


Proof. Since the norm || • || is induced by an inner product, it is smooth away 
from 0 (cf. the discussion in [27, Proposition 13.14]). In particular, this entails 
that the scaling maps (5) are smooth by [27, Theorem 42.13]. 

To see that 1Z is a diffeomorphism, observe that T: C°° (/, g) —>■ C™ (/, G) is a 
(smooth) diffeomorphism by [17, Theorem A]. Its inverse is the right logarithmic 
derivative 6 r : Cf°(I,G) -» C°°(I,g). Hence 7Z = sc -1 o 5 r | Imm (/,G) n c“(J,G) is 
a diffeomorphism whose inverse clearly is 7Z~ X . As 7 Z~ l = To sc is also smooth, 
the assertion follows. □ 

The right action V x Diff + (/) —> V, (c, ip) i —> c o <p restricts to an action 
on T 5 *. We can now define a reparametrization invariant (pseudo)metric 5 on V * 
as follows. 

Definition 3.10. Restrict d-p to a pseudometric dp t on V*, i.e. from (4) we 
derive 

dvA c o,ci) ■■= j |ko(i) - qi{t)\\ 2 dt = d L 2(1Z(c 0 ),1Z{ci)), (7) 

where qi := 7 Z(ci), i = 0,1. 


3.4 Riemannian geometry of the SRV transform 

There is a geometric interpretation of the distance function dp t obtained via the 
SRVT that motivates our choice of approach. In the present section we explore 
this interpretation in the context of Riemannian geometry on spaces of curves. 
Consider the space of curves C°°{I , g) with the L 2 -metric 


( h,g ) L 2 := J{h(t),g(t))dt. 


This metric defines a (weak) Riemannian metric 6 on C 00 (/, g). Moreover, since 
the image of the SRVT is an open subset of C°°(/, 0 ), the Riemannian metric 


5 We will see later see Theorem 3.11 and Theorem 3.16 that the pseudometric is the geodesic 
distance with respect to a first order Sobolev metric if dim g > 2. Hence in these case we 
actually obtain a metric on V *. 

6 Here the term “weak” means that the Riemannian metric does not determine the topology 
on T c C°° (I, g) = C°° (7,g). 
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restricts to a Riemannian metric on the image of the SRVT. Now we exploit 
that the SRVT is a diffeomorphism (Lemma 3.9), to obtain a pullback metric 
on TV whose geodesic distance will turn out to be d-p, ■ Before we prove this, let 
us derive a formula for the pullback metric. 

Recall that S r c := R~»(c) is the right logarithmic derivative whose tangent 
map at c we denote by T c S r . 

Theorem 3.11. Let c £ V* and consider v,w £ T C V *, i.e. v,w: I —> TG are 
curves with v(t),w(t) £ T C ^G. The pullback of the L 2 -metric on C°°(I, g\{0}) 
under the SRVT to the manifold of immersions V * is given by: 

G c (v,w)= f j (D s v, u c ) (D s w, u c ) 

J i (8) 

+ ( D s v - u c ( D s v , u c ) , D s w - u c ( D s w, u c )) ds, 

where D s v := T c i5 r (u)/||c|| ) u c := <5 T ’(c)/||<5 r '(c) || is the unit tangent vector of 
S r (c) and ds = ||c(f)|| dt. Consequently , the pullback of the Lr-norm is given by 


f 1 2 2 

G c (v,v) = J -(D s v,u c ) + \\D s v - u c (D s v,u c )\\ ds. 


Proof. We have to compute the tangent map of the SRVT at a curve c in the 
direction of a vector field v along c. Recall that 1Z = sc o <5 r , whence the 
chain rule implies T c lZ(v) = T sc o T c S r (v). Hence setting z v := T c 5 r (v) £ 
Tr< c) C°°{I, g \ { 0 }), it suffices to compute TV^sc^.y). 

By Proposition 3.2, there is a smooth map 7 : ] — e, e[xl —t 0 \ {0} with 
7 ( 0 , t) = 5 r (c)(t) and ^7 (x,t) = z v (t). Thus the tangent map of sc can be 
computed as follows. 


T c TZ(v) = 7V( c )Sc(Zy)(f) = 


dx 


sc(j(x,t)) = 


x—0 


dx 


li x i t) 

c=o V\hi x ,t)\\ 


= \\6 r c(t))\\-h v (t) - -rc(t))||-i (z v (t),6 r c(t))S r c(t). 


Set z w := T c 5 r (w) and substitute into the definition of the pullback metric these 
formulae. Then we simplify the expression as follows: 

G c {v, w) = (T c TZ(v),T c TZ(w)) l2 

= j rcWII- 1 (z v {t),z w {t)) - |||5 T 'c(t)||- 3 < z v {t) 1 8 r c(t)) (z w (t),S r c(t))dt. 

(9) 

We set u c {t) := ||^ r ^*j|| and notice that ||it c (f)|| = 1 . Inserting these identities 


12 







in (9) we can simplify G c (v,w) as follows 


J n^)ir i «^( i ),^(i)) - 4 rc(i)r 2 <*,<*), rcw) <z w w,^cw))dt 
= J IKc^ir 1 ^z v (t),z w {t)-^u c (t) (z w (t),u c (t))^ dt 

= j i ll^c(i)||^ 1 Q< z v (t),u c {t))(z w (t),u c {t))+ 

+ (z v (t) - u c (t){z v (t),u c (t)),z w (t) - u c (t)(z w (t),u c (t)))^ dt 
f ncr , sn ( 1 / T c S r (v)(t) , A /T c 6 r (w)(t) . A 

“/, rc(t)l1 (i (■W tt ' (t) ) (wF’ ^V + 

/ ^ / z v (t) n\\ z w(t ) ^ / z w (t) 


11^)11 


11^)11 


11^)11 


\\S r c(t)\l 


dt 


In passing from the second to the third line we have used the trivial identity 

z w (t) = u c (t)(z w (t),u c (t)) + (z w (t) - u c (t)(z w (t),u c (t))). 

We will now exploit that the norm and the Riemannian metric on the tangent 
spaces are right invariant, i.e. they are invariant under right translation. This 
entails ||(5 I '(c)|| = ||c||, whence ds = ||c(t)|| d£ = 11<5 T ’c(t)11 dt. Substitute this 
together with the identities for z v ,z w and u c to obtain 

G c (v,w) = J ^(D s v,u c ) (D s w, u c ) 

+ ( D s v - u c ( D s v , u c ), D s w - u c ( D s w , u c )) ds. 

Thus (8) holds and the formula for the norm follows directly by specialization. 

□ 

Remark 3.12. For computations in the rest of the paper, we will need neither 
the explicit form of the pullback metric on V * nor the norm computed in Theo¬ 
rem 3.11. The idea is to use the SRVT to relegate all questions concerning the 
metric to (C°°{I, g \ {0}), (•, ■)&). 

Remark 3.13. Theorem 3.11 generalises [3, Theorem 4.2] (for the parameter 
values a = 1 = 2b): Viewing the vector space R 2 as an abelian Lie group, our 
approach recovers the constructions of the pullback metric in the vector space 
case. 

Notice however that we take a slightly different (but equivalent) perspective 
on the shape spaces involved in the construction: In [3] the unparametrized 
curves are modelled as the quotient of Imm(/,1 2 ) modulo translations. Pick¬ 
ing representatives for each class, one can show that on the level of infinite¬ 
dimensional manifolds this yields the same concept of (unparametrized) curves 
as our approach. Moreover, as the SRVT is translation invariant (Lemma 3.6), 
the two constructions yield the same Riemannian manifold (as already observed 
in [3]). 
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The pullback Riemannian metric constructed in Theorem 3.11 defines a Rie- 
mannian metric on the space of parametrized curves V *. For curves which take 
their value in R™ the two terms in the integral can be seen as measuring bending 
and stretching deformations, respectively [3, 42]. Therefore, metrics of the form 
(8) are known as elastic metrics. 

Note that the formula for the Riemannian metric is given in terms of a 
tangent map of the right-logarithmic derivative. We here give an explicit formula 
for this map. A proof can be found in Appendix A.l. 

Proposition 3.14. Let c £ C°°(I,G) and v £ T c C°°(I,G), i.e. v £ C°°(I,TG) 
with v(t) £ T c ( t )G, Vt £ I. Then writing [•, •] for the Lie bracket in g we have 


T cS r (v)(t) = («(*))) 




In what follows, we discuss the Riemannian geometry of the L 2 -metric on 
the image of the SRVT. This will shed light on the geometry induced by the 
pullback metric on 'P*. Geodesics in the image of 7Z are just restrictions of 
geodesics in g) (with respect to the L 2 -metric) to Im 1Z = C°°{I, g\ {0}). 

These geodesics correspond to geodesics in V * with respect to the pullback 
metric (8). As we are interested in the geodesic distances on V*, let us first try 
to understand geodesics on ImTJC g). 

Proposition 3.15. Consider C°°(I, g) with the weak Riemannian structure 
induced by the L 2 -inner product ( f,g)L 2 : = h(f(t),9(t))dt. 

1. The space (C°°(I, g), flat the sense of Riemannian geometry. 

2. The open subset C°°{I, g \ {0}) C C°°(I, g) is also flat. Furthermore, 

for dim g > 1 there exist points in g \ {0}) not connected by a 

minimizing geodesic. 


Proof. 1. The Riemannian metric is the inner product of a vector space. 
Hence g) is flat as all derivatives of the Riemannian metric with 

respect to the base point vanish. 

2. Note that C°°(I, g \ {0}) is an open subset of C°°(/,g), whence flat. For 
dim g > 1 the set g \ {0} is connected but not convex, whence there are 
points not connected by minimizing geodesics, e.g. for c £ g \ {0}) 

the minimizing geodesic connecting c and —c is not contained in C°°(I, g\ 
{0}) (since c(t) + |(— c(t) - c(t)) = |c(f) + |(-c(t)) = 0). □ 

As the Riemannian structure of C°°(I, g\{0}) is induced by C°°(I , g), the L 2 - 
distance locally describes the geodesic distance. Thus the distance function (8) 
locally describes the geodesic distance with respect to the pullback Riemannian 
metric on V *. Since minimizing geodesics need not exist between points f,g£ 
C°°(I,q \ {0}), the geodesic distance might be strictly greater than the L 2 - 
distance if f(s) + t(f(s) — g(s )) = 0. 

However, we will prove that, at least if the dimension of the Lie algebra is 
large enough, the L 2 -distance coincides with the geodesic distance (see Theorem 
3.16). To prove that the geodesic distance on G°°(J, g \ {0}) coincides with the 
L 2 -distance, we have to approximate the minimizing geodesic by paths in g\{0}. 
To avoid a lengthy exposition at this time, we have relegated these details to 
Appendix A.2. 
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Theorem 3.16. If dim g > 2 then the geodesic distance of C°°(I, g \ {0}) is 
globally given by the L 2 -distance. In particular, in this case the geodesic distance 
of the pullback metric (8) on P* is given by the distance function (7). 

Remark 3.17. Notice that Theorem 3.16 entails that the pseudometric dp* 
defined via (7) is non-degenerate for dim g > 2, i.e. d-p* (c, c) > 0 if c ± c.. Thus 
dp* is a metric on P*. In particular, we will say that dp* defines a distance 
function on P* in this case. 


The crucial observation here is that one obtains a distance on P*. This 
enables us in Section 5 to compute distances between curves and to deform 
curves into each other along geodesic paths. The relation of geodesic distance 
and L 2 -distance has previously been used in animation classification tasks in 
for d > 2 e.g. in [15]. 

Finally, we can follow the argument given in [3, Theorem 6.1] to derive 
information on the curvature of the shape space: 

Corollary 3.18. The curvature of the space 5* := P*/Diff + (J) with the Rie- 
mannian metric induced by the pullback metric (8) is non-negative. 


Proof. The Riemannian structure on 5* is induced by the one on P*, i.e. the 
canonical quotient map P* —> 5* = P*/Diff + (/) is a Riemannian submersion. 
Hence we can apply the O’Neil curvature formula (see e.g. [11, Theorem 3.20]) 
for the sectional curvature of the quotient 5*. For orthonormal vector fields 
X , Y on S* this yields: 


Ks.(X.Y) 


Kp*(X,Y) + ± 


vert 


x,y 


2 


Here X and Y are horizontal lifts of X and Y to T* and 


vert 


x,y 


denotes 


the vertical projection. As V* is flat by Proposition 3.15, the curvature Kp* 
vanishes, whence I\s* is non-negative. □ 


3.5 Distance as an optimization problem 

In the last section we acquired a distance function on P* for dim g > 2. Using 
this distance function, we can calculate distances in the shape space := 
P*/Diff + (/) 7 by solving the optimization problem 

ds. ([<*], [ci]) := inf ( [ \\q 0 (t) - gi(v?(<))\/^(t)|| 2 dd , (10) 

veDiff+p) Vd/ J 

where qi := TZ(ci), i = 0,1. 

When computing d$* we therefore need to perform an optimization over the 
diffeomorphism group Diff + (I) (see for example [31] for more details and an 
extension to piecewise linear curves instead of immersions). 

In practice, one of two different algorithms is used to solve this optimization 
problem: either a gradient descent based approach or a dynamic programming 

' Recall that in the beginning of the Section 3 we set out to construct a distance function for 
the shape space S = 'P/DifUf /). As explained in the last section, this is not possible within 
the SRVT-framework. Hence we have to use the smaller space 5* (which can be identified 
with a subset of 5). 
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(DP) algorithm. For our numerical experiments and applications in Section 
5, we have used DP, which constructs a piecewise linear approximation of the 
optimal reparametrization p in (10). See [6, 38] and references therein for more 
information on the use of DP for shape analysis. 

4 Closed curves 

In various applications, we are particularly interested in closed curves. For 
example, in object recognition, closed planar curves can be used to represent 
outlines of objects. In computer animation, closed curves are cyclic animations 
that can be repeated multiple times with no visually noticeable discontinuities. 

Here we will derive a method to calculate a closed curve approximation 
of an existing open curve. Again, we will be working with SRV transformed 
representatives q of curves c G V *. We denote the image sets via the 1Z transform 
of the open and closed immersions on G by 

C° := K(lmm(I, G)) = g \ {0}) = {q G C°°(I, g) \ q(t) ^ 0 Vi £ 1} 

C c := 7 Z({c € Imm(/, G) \ c( 0) = c(l) = e}), 

respectively, i.e., we will be working exclusively with parametrized curves, and 
not with shapes. Notice that the open curves are an open subset of the locally 
convex space C°°(I,g). 

A curve c in the set of closed curves C c will in general only be C°-closed. 
This means, we do not require that the derivatives ^c(O) and ^c(l) (or higher 
derivatives at the closing points) coincide. Hence the closed curves considered 
here admit “corners” at the closing points and should not be confused with 
smooth loops in G , i.e. smooth maps from the unit circle S 1 to G. 

Consequently, the closed curves computed via the methods in this section 
will only be C 0 -closed. Though a higher order closing might be desirable, for 
the applications we have in mind the closing of the curves is sufficient. 

4.1 Closed curves as a closed submanifold 

We will now prove that C c is a closed submanifold of the open curves C° and 
thus a closed submanifold of C°°(I , g). Since 7Z is a diffeomorphism, this entails 
that {c e Imm(I, G) | c(0) = c(l) = ej = 7\1 _1 (C C ) is a closed submanifold of 
P* = {/eImm(/,G) |/(0)=e}. 

The basic idea is now to construct C c as the preimage of a closed submanifold 
under a submersion. To this end, consider the point evaluation map 

ev x : C°°(/, G) —> G, f /(1), (evaluation in 1). 

It can be observed that 

C c = r^ 1 (e), r: C° —> G, r := evi o T o sc, 

where sc denotes the scaling map from (5) and we note that r(q) = 7?. _1 (g)(l). 

Proposition 4.1. The map r: C° —» G is a submersion. Hence C c = r -1 (e) is 
a closed submanifold of finite codimension in C° C C°°(I,g). 
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Proof. Let us first establish the smoothness of r. In the proof of Lemma 3.9 we 
have already remarked that L and sc are (smooth) diffeomorphisms. Further, 
evi is smooth as a consequence of [27, Theorem 42.13] whence it restricts to 
a smooth map on the closed submanifold C£° (/, G). By abuse of notation we 
will denote the induced map on the submanifold also by evi. We conclude that 
r is smooth and in particular it will be a submersion if and only if evi is a 
submersion. 

For general locally convex (infinite-dimensional) manifolds a map is called 
submersion (see [18]) if for every point in its domain we can find submersion 
charts, i.e. charts around the point and its image which turn the map locally 
into a projection . 8 We will now construct submersion charts for c £ C%°(I, G). 
To this end, let ip '■ G A V —>• W C g be a chart for G around the identity. Then 

ip *: c ?(/, G)nc°°(/, v) -> {h e c°°(i , fl ) | h( o) = o}nc°°(/, w), 

is an identity chart for G). Since C*°(7, G) is a Lie group, translating ip* 

by c yields a chart centered at c. Denote by evi: {h £ C°°(I, g) | h{ 0) = 0} —>■ 
g, h i->- h( 1) the point evaluation in 1. We compute for v £ {h £ g) | 

h(0 ) = 0} 

evi o L c o ip^pv) = evi o L c (i/f _1 o v) = c(l)^ _ 1 (u(l)) = (L c(1 ) o ip~ 1 )ev 1 (v). 

As L c ,L c t i) and the charts if, ip* are diffeomorphisms, it suffices to construct 
submersion charts for evi. To construct these charts, we split {h £ C°°(I, g) | 
h{ 0 ) = 0 } non canonically: 

A: {h £ C°°(/,g) | h(0) = 0} -> {h £ G°°(J,g) | h{ 0) = h{ 1) = 0} x g, 

h ha (t h(t) — th( 1 ), h( 1 )) 

Since evi is smooth, A is smooth. Furthermore, A is a isomorphism of locally 
convex spaces as A _ 1 (c, q) = (f £ c(t) + tq ) is also smooth. Now let pr 0 be 
the projection onto g in the above product. Then one computes evi o A = pr g . 
Hence, A and id 0 form a pair of submersion charts for evi, whence evi and evi 
are submersions. We deduce that r is a submersion. 

To prove the final assertion we use that C c = r^ x (e) is the preimage of a 
point under the submersion r. Invoking a version of the regular value theorem 
[18, Theorem D] for infinite-dimensional manifolds, we deduce that C c is a closed 
submanifold of C° (whence also of C'°°(/,g)). If G is finite-dimensional, C c is a 
submanifold of finite codimension. □ 

Remark 4.2. The right Lie group action V * x Diff + (/) —» V*, (c, ip) >—»■ c o ip 

restricts to an action on the closed curves {c £ Imm(/, G) \ c(0) = e = c(l)}. 
Hence the pullback Riemannian metric ( 8 ) induces a reparametrisation invariant 
Riemannian metric on the space of closed curves (equivalently on C c ). We 
remark that C c will in general not be flat in the sense of Riemannian geometry 
(cf. the computation of curvature in [3, 5.2]). 

Remark 4.3 (The framework for infinite-dimensional Lie groups). If we assume 
that G is a Hilbert Lie group, the results obtained so far carry over without any 

8 In the case of a finite-dimensional Lie group G, [18, Theorem A] asserts that a sufficient 
condition for evi : GJ° ( I. G) —> G to be a submersion is that its differential at every point is 
surjective. This is easily verified. 
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changes in the proofs. Here by Hilbert Lie group we mean a Lie group in the 
sense of [27, Section 36] or equivalently [35] modelled on a Hilbert space. 

More generally, some of the results are still valid for Banach Lie groups 
which are modeled on a Banach space with a smooth norm away from 0 (see 
[27, Chapter 13] for more information on smooth norms). Notice however that 
in this case complications arise. For example, in the Banach Lie group case we 
lack an inner product. Hence, one has to replace the formula for the derivative 
of the norm in the proof of Theorem 3.11 and it is unclear how to relate the 
distance d-p t to the geodesic distance from Riemannian geometry. 

We address now the problem of projecting open curves of C° onto the sub¬ 
manifold of closed curves C c . In theory, we could define a projection from C° 
onto C c by stating a constrained minimization problem 



where qo is the curve to be approximated. Instead of minimizing the distance 
from closed curves to qo in what follows we opt for minimizing the closure 
constraint. This approach leads to the formulation of a gradient flow whose 
solution is the desired closed curve. 

4.2 Projection via a gradient flow 

Given an SRV representative q of an open curve in the Lie group G, i.e., q £ 
C° \C C , we will try to “close” the curve by enforcing (at least approximately) 
the closedness constraint r(q) = TZ~ 1 (q)(l) = e. 

We can measure the distance between the identity and the actual endpoint 
by using the functional 



( 11 ) 


where log denotes the inverse of the exponential map exp : g —> G, and is defined 
in a neighbourhood of the identity of G. Notice that 



Note also, that this only measures first order continuity. 

In many practical applications, we are satisfied with $>(q) ~ 0; particularly 
in the computer animation applications discussed later in Section 5 a precise 
enforcement of the constraint < f>(g) = 0 is not necessary to achieve visually 
pleasing results. This leads us to the idea of closing curves by minimizing the 
functional $, and since C° is an open subset of a vector space, we can use a 
straightforward gradient descent method for solving this problem. 

In order to compute the gradient of 4>, we will need both the tangent map 
of the evolution operator T, for which we get [27, 38.10 Corollary]: 



( 12 ) 


as well as the following tangent map for TZ 1 



(13) 
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Here, Ad g denotes the adjoint representation of the Lie group: Ad g = L gir o R~f. 
Since T q r = T q ev i o 7W 1 the formula (13) immediately yields a formula for 
T q r. Moreover, as the tangent map of evi can be canonically identified with an 
evaluation in 1 (cf. proof of Proposition 4.1), we see that T q r is given integrating 
to t = 1 in the formula (13). 

Furthermore, we need the tangent of the logarithm map. To this end, let us 
first consider the tangent of the Lie group exponential map T exp: Tg = g x g —► 
TG. Using the right trivialisation of the tangent Lie group TG = g XAd G, we 
have 

T u exp(u) = J f ? exp( „ ) * o dexpju), (14) 

for a unique (linear) map dexp u : g — > g called the right trivialized tangent of 
the exponential map. Following [9, 20] we remark that dexp M satisfies: 


dexpju) 


n+ ^\u,v\ + i[u, [u,u]] + ... 


oo 


0 )> 

z—&d u 


(15) 


where ad M is the derived representation of the Lie algebra ad u v := [w, v]. By [35, 
Remark II.5.8] the formula (15) holds also for infinite-dimensional Lie groups 
modelled on Banach spaces. Note that dexp M (u) = u and dexp u = 5 r (exp) u (cf. 
[35, p.340 - 341]). 

By using that expo log = idc in a neighbourhood of the identity in G, (14) 
yields 


T g log = dexp log , -j (R g - 1*) , for all g in the domain of log. (16) 

Combining these tangent maps, we can calculate the gradient of the error func¬ 
tional <F in (11). We use the following notation for adjoints of bounded linear 
operators. 

Definition 4.4. Let A: g —>• g be a bounded linear operator and (•,•) be an 
inner product on g. Then we denote by A 1 the adjoint operator of A, i.e. for all 
v,w £ g we have (A(v),w) = (v, A'(w)). 

Theorem 4.5. The gradient of the the error functional 

<F : C° —>■ R, $(g) := i|| log (^~ 1 (g)(l)) || 2 , 
with respect to the L 2 inner product is the vector field on C° given by 

grad(4>)(g) = ||g|| 9(q) + (9(q), Tpy) q, (17) 


where 


9(q) := Ad^ Ad* (?) (dexp log(r(9)) ) 1 ’ (log(r(g))) £ G°°(J,g) 
c(q) ■= R~\q) £ C?(I,G) and r{q) := ^ 1 ( 9 )( 1) £ G. 
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Proof. Let / be a vector field along the curve q. The L 2 gradient of $ is then 
defined by 


T q ®(f) = (grad($)(<?),/) l 2 = J (grad(<f>)(g), /) dx, 


and 


T q *(f) ■■= 


de 


e— 0 


log^-^g + c/)(l)), log^fo + e/)(l)) 


By differentiation we obtain: 


T q $(f) = (T q (logor)(f), (logo r )(q)) 

( = ) ( dex Piog(r( q )) Rr(q)-uT q r (/), log(r(g))) 

= (T q r (f),Rt (q) -iM ex Piol(r( q )))' l °S(r(q))) ■ 


Insert now the formula for the tangent maps (12) and (13) and observe that the 
integral appearing is a weak integral, i.e. we may interchanged the integral with 
any of the continuous linear functionals (-,v). Hence the above formula for T q ' I> 
simplifies as follows: 


T q®(f) = ^L r ( 9 )*^Ad c(5)(l) -i ^/(a:)||g(a;)|| + (f(x), q(x)) dx > 

R l( q )~ 1 * ( dex Piog(r( 9 )) ) f log(r (q )) ^ 

= ^ ^Ad c( , )(x) -i (j(x)\\q(x)\\ + {f(x),q(x)) > 

L t( q )* R i(q)-i* ( dex Pk>g(r( g ))) f lo g(r(<?))^ dx 
= ||g(a:)|| + {f(x),q{x)) ^ ,0(q)(x)j dx 

= j^f(x)\\q(x)\\,9(q,x)) + (f(x),q(x)) (^ ^J^ ,9(q)(x)^ dx 

= j i (\fix), 0(q)(x)\\q(x)\\ + (e(q){x), Q( x Y) dx » 

Reading off the gradient of $(g) from this expression we derive: 

(gra,d($)(q),f) L2 = (e(q)\\q\\ +(6{q),^)qJ^ . □ 

The projection onto the space of closed curves C is then obtained by solving 
the differential equation for u(t,r ): 

du 

— = -grad($)(u), u(t, 0 ) = q(t). 
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Remark 4.6. The results of the present section on tangent maps carry over 
verbatim to the case of infinite-dimensional Lie groups modelled on Banach 
spaces. This is due to the fact that Banach Lie groups are locally exponential 
Lie groups (i.e. the Lie group exponential map is a local diffeomorphism near 
the unit) cf. [35, Proposition IV.1.2]. Moreover, the formula for the trivialised 
tangent map of the exponential still holds in this case by [35, Remark II.5.8]. In 
particular, this shows that Theorem 4.5 remains valid for Lie groups modelled 
on a Hilbert space. 

The derivative dexp u of the Lie group exponential appears in the formula 
(17). If we want to apply the formula for the gradient in computations we thus 
need to compute this derivative. However, for certain finite dimensional Lie 
groups, these additional computations can be avoided as the derivative vanishes 
in the formula. We will prove now that in 

%) = Ad I( g )-i Ad r( 9 ) ( deX Plog(r( 9 ))) t MK?))) 

the term (dexp[)) g , ( ? p)^ vanishes if we choose the inner product induced by the 
Cartan-Killing form of a compact and semisimple Lie group. 

Definition 4.7. Let G be a finite-dimensional Lie group with Lie algebra g. 
Define the Cartan-Killing form 

gxg->R, Kg(x,y) ■= tr(ad x o ad y ), 

where we denote by tr the trace of a linear map. Recall from [21, Remark 
12.2.14] that n g is a negative definite form if and only if G is a compact and 
semisimple Lie group. 

Corollary 4.8. Let G be a compact and semisimple Lie group. With respect to 
the L 2 -metric constructed from the inner product —n g induced by the Cartan- 
Killing form, the gradient of the the error functional $ (11) is the vector field 
on C° given by 

( 18 ) 

where 

a{q) := Ad^-i Ad* (g) (log(r(<?))) £ C°°(I,g) 
c{q):=n~ 1 {q)€C^{I 1 G) and r(q) := TZ~ 1 (q){l) £ G. 


grad(<f>)((?) = HgH a(q) 


Proof. As a shorthand we define u := log(?'(g)). Having the formula (17) at 
our disposal, it suffices to prove that the adjoint operator (dexp,” 1 )'!' fixes u. 
As (dexp“ 1 )i' = (dexp[j) -1 it suffices to prove that dexp^ fixes u. Recall from 
[21, Exercise 5.4.5] that the operators ad x are skew symmetric with respect to 
the Cartan-Killing form, i.e. K g (ad s (y ), z) — n g (y, —ad x (z)) or in other words 
(ad x )^ = — ad x . Using the series identity (15) for dexp„ together with the fact 
that the mapping A is a continuous algebra morphism, we obtain the 

identity 


/ OO 1 \ oo 

(dexpj t = ( E (fcTl)! (ad “ )fc I = £ 


i /c=0 


°° 1 

(adt) fe =E77^T(-a d .) fe 


fc =0 


(fc + 1)! 


k =0 


(k + 1)! 


Hence ad u (w) = 0 implies (dexp u p(u) = u and the assertion follows. □ 
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5 Numerical Results 


In this section we present some results obtained with the methods developed in 
the previous sections. 

First, in Section 5.2, a few simple examples of curves in 50(3) will serve 
to demonstrate the basic ideas discussed in this paper: interpolation between 
curves and closing of open curves. Next, we will present some numerical results 
for the specific applications to computer animation problems: interpolation be¬ 
tween existing motions in Section 5.3 and removing discontinuities in (almost) 
periodic motions in Section 5.4. 

5.1 Some implementation notes 

The motion capturing data used in the two animation application Sections has 
been taken from the CMU motion capture database [8]. 

Computational tools for 50(3). In the examples we consider the Lie group 
50(3) (or a product of multiple copies of 50(3)), for which we will use a matrix 
representation, so left and right translations correspond to matrix products. 
The Lie algebra so (3) consists of the 3x3 skew-symmetric matrices which are 
isomorphic to vectors in M 3 via the hat map: 


1 X1 1 

k j 

f 0 

-xa 

X 2 \ 




0 

—X\ 

V ^'3 ) 

' ' 

{ -X 2 

Xl 

o / 


Efficient ways to compute the exponential map and the logarithm in 50(3) 
are available in the literature, see [9, 10, 22], In this setting, the Lie algebra 
exponential exp : so(3) —>• 50(3) can be efficiently computed using Rodrigues’ 
formula: 


exp(i) = I + 


sin(0). , 0 sin 2 (0/2), 2 

e + 0 2 ’ 


0:=M- 


Similarly, there exist efficient means to compute the logarithm of an orthogonal 
matrix X G 50(3): 


l°gPO 



X ^ /, and X close to /, 


where y = |(X — X^). 

As in the previous section we denote here by A'f the adjoint linear operator 
(i.e. the transpose of the matrix A'). 

The Lie group 50(3) is compact and semisimple (see [21, Lemma 2.1.4 and 
Example 5.5.4]). To compute the gradient of the error functional 4) (11) as in 
Corollary 4.8, we will thus assume that the inner product on g is the negative of 
the Cartan-Killing form n B . For 50(3) it is well known that the Cartan-Killing 
form is given by 

K„( 3 ){X,Y):=tr(XY) 

As X, Y are skew symmetric, the corresponding inner product turns out to be 
the familiar Frobenius inner product — k S o( 3 )(X, Y) = tr(A'Fl). Notice, that 
one can proceed similarly if G is a (finite) product of copies of 50(3). 
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Discrete curves. Given a continuous curve c in SO( 3), we approximate it by 
picking a discretization {Oi}™ =0 of I and constructing a curve c based on discrete 
points { Ci := c(#i)}" =0 , between which we interpolate along geodesics in SO(3): 


71—1 


C W : = VX[0 fc A + i)(O exp 


k =0 


£ — i 

h+1 


— Oi 


■log(c fc+ ic[) ) Cfc, 


(19) 


where x is the characteristic function. 

Computing the square root velocity transform (3) of such discrete curves 
then results in piecewise constant functions q = {qi}^!^ in the Lie algebra, with 
discrete points qt given by: 


- Vi ._ log(ci + icj), 

dH’ »*«-»* ' 

The inverse SRVT (6) of a piecewise constant function q in g, is given by a 
piecewise geodesic curve c in G, as formulated in (19), with points c, given by: 

c i+ i = expfl|gj|| qi)ci, i = l,...,n-l, c 0 = e. 

Curve reparametrization. As we have seen, reparametrizations in Diff + (J) 
act on curves from the right. For discrete curves, this means either a change 
in the underlying grid or a change in the discrete points. For our numerical 
experiments, we have chosen to keep a fixed grid and resample curve points. 
Applying a reparametrization ip £ Diff + (/) to the discrete curve c then results 
in a new discrete curve c with sampling points {ci}" =0 computed using geodesic 
interpolation: 

ci := Cj-exp(slog(cj + ict)), s := — p-, i = 0,...,n, 

where j is an index such that 9j < p(0i) < 9j+ 1- Note that cq = cq and c n = c n 
naturally follow from the definition of the diffeomorphism group Diff + (/). 

Curve interpolation. To perform interpolation between two parametrized 
curves Cq and ci, we interpolate linearly between their SRVT representations 
and reintegrate the result: 

[ 0,1 ]xP.xP.-K%°(I,G) 

(s, Co, Cl) I-X ^ _1 ((1 - s)ft(co) + 371(a)), 

with interpolation parameter s. The construction of the curve is possible due 
to the vector space structure of the Lie algebra g and yields precisely a geodesic 
of G°°(J, g) (cf. Proposition 3.15). Notice that the interpolation map takes its 
image in C£°(I,G) and not in V* , i.e. in general the result will not again be 
an immersion. This is due to the fact that elements in V * take their image in 
g \ {0} and this set is not convex in g. In practice this problem is not very 
serious, however it seems difficult to exclude this problem without turning to a 
cumbersome set of conditions on the initial data. 

Curve closing. As we have already mentioned the Lie groups in our main 
example are compact and semisimple, whence Corollary 4.8 is applicable. In 
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particular, the gradient for the curve-closing method in Equation (18), takes 
the following form for curves in SO( 3): 

grad(<f>)(g) = ||g||clog(c(l))c t + (clog(c(l))c t , r^) q 

IMI 

where c := TZ^ 1 (q) (computed via Lie-Euler integration), and the gradient flow 
can be discretized as 


u k+1 =u k — Q!fe grad($)(u fc ), 

where every u k is a discrete curve as defined above, i.e., u k = {u.f}" = 0 . Note 
that for the curve closing, a curve’s parametrization is fixed. 

This iterative approach allows us to balance accuracy and computational 
expense, which is useful in the computer animation applications discussed in 
Section 5.4. 

5.2 Curves on SO{3) 

We start with two simple demonstrations of the methods developed in this 
article: one interpolation between curves on SO( 3) and one application of the 
curve-closing algorithm laid out in Section 4. In order to visualize a curve 
c : I — > SO(3), we will take one or more unit vectors in M 3 , rq,..., v n , and plot 
the curves resulting from transforming these vectors by the successive elements 
in SO( 3), i.e., we will plot Vi(t) := c(t)vi for i = 1,..., v n . The curves fq(t) will 
therefore evolve on the unit sphere in R 3 . 

Figure 2 shows the result of interpolating between curves ci, c% : I —> SO( 3) 
using the approach outlined in Sections 3 and 5.1. 

This is the basic mechanism underlying the shape distance computation: 
We represent shapes, i.e., equivalence classes of parametrized curves under 
reparametrizations, by a representative, i.e., a single parametrized curve. Then, 
given two shapes and two corresponding parametrized curves, we try to find a 
parametrization of one of the curves that minimizes the distance to the other 
curve. This minimum distance is then, according to Lemma 3.4, the distance 
between the two shapes, i.e., the equivalence classes. 

The two curves in the first row of Figure 2 represent the original curves 
between which we interpolate using Equation (20). In the middle figure, we 
see their parametrizations - whereas the orange dashed curve has a uniform 
parametrization, the blue one is in a sense compressed in the beginning and 
then stretches out. When interpolating between those two parametrized curves, 
as seen in the second row, the resulting interpolation first contracts and then 
expands. In the right side figure in the first row, the lower curve has been 
reparametrized, using dynamic programming, to better match the two curves. 
This minimizes the distance between the two curves as discussed in the previous 
section, and the third row of figures shows the corresponding interpolating path. 


In Figure 3, we see the results of applying the curve closing algorithm to an open 
curve in SO( 3). We have plotted the results at different stages of the iterative 
algorithm to highlight the evolution of the closing process, i.e., how the curves 
iteratively move towards closedness. In this toy-example, we have chosen a 
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Figure 2: Interpolation between two curves in SO( 3) with and without 

reparametrization. The curves shown here are the points traced out by the 
vector (1, 0,0) T in R 3 when transformed along the curves in the Lie group. The 
top row shows the original curves, their parametrizations and a reparametriza¬ 
tion performed to minimize the distance between the curves. The thick red 
curves in the second and third rows are points along the geodesic path between 
the two original curves, at times s £ {1/4,1/2,3/4}, from left to right. 
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Figure 3: Application of the closing algorithm to a single curve in SO( 3). The 
two curves shown here are the points traced out by the vectors (1, 0, 0) T and 
(0,1,0) T in M 3 when transformed along the curve in the Lie group. The top 
left figure shows the open starting curve, whereas the remaining panels show 
the evolution of the closing algorithm. In the last panel, the original curves 
are superimposed as dashed lines to show the deviations caused by the closing 
method. Note that for visualization purposes, a small stepsize was chosen, 
resulting in more iterations than would otherwise be necessary to achieve a 
satisfying accuracy. 
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Figure 4: Interpolation between two walking animations, shown in the top 
row. The two lines shown in the plots are the trajectories of both feet, with the 
big steps over the obstacles clearly visible. Below, interpolations between the 
two animations computed using different methods are plotted. The upper row 
shows existing results from [6], where additional feature point information was 
used to guide the interpolation. The bottom row shows interpolation results for 
matching the SO(3) d animations without (on the left) and with (on the right) 
reparametrization. No additional feature point information was necessary for 
this method. (Note that the feature matching approach by design requires 
reparametrization.) 


small stepsize to better show the behaviour of the algorithm. In practice, just 
two to three iteration steps are typically enough to achieve a sufficient degree 
of closedness in the curve. 

5.3 Motion interpolation 

We now want to apply the interpolation method outlined earlier and demon¬ 
strated in the previous Section to entire motions of virtual characters. This 
means that instead of curves in SO( 3), we are now dealing with curves in 
SO( 3) d , with a copy of SO( 3) for every joint in the animated character. The 
numerical methods as described in Section 5.1 stay the same however. 

In Figure 4 we start with two original motions: One with the animated char¬ 
acter stepping over a high obstacle and one with the character stepping over a 
low obstacle. We want to interpolate between those two motions, allowing us 
to let the character step over obstacles of arbitrary heights. This problem was 
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already considered in [6] , where it was found that the approach of parametrizing 
animations with Euler angles and then performing shape analysis [15] produced 
unsatisfactory results for these motions and had to be augmented with extra 
landmark information to achieve a realistic interpolation result. Here we find 
that, for this problem, the use of a Lie group formulation allows us to achieve 
similarly good results as in [6], but without the need for additional landmark 
information. Figure 4 shows results for the Lie group formulation once without 
and once with reparametrizations. Reparametrizations in the context of anima¬ 
tions, also known as time warps [24], serve to align two animations on the time 
axis, i.e., they can speed up or slow down parts of an animation to match the 
motions more closely. We see that in this case, reparametrizations are necessary 
to achieve satisfying blending results. This is likely due to a slight phase shift in 
the beginning (compare the starting points of the trajectories on both original 
animations in Figure 4). 

A natural extension of this scheme would be to interpolate between multiple 
animations. Among others, this could be useful for extended walking animations 
where combining multiple animations and varying the interpolation weights over 
time could help produce a large but consistent set of motions to avoid repetitive 
visual elements. This could for example be accomplished by formulating this 
animation blending problem as a Karcher mean ([42, 43]) of multiple animation 
curves. We will consider such a scheme in future work. 

5.4 Periodic motions 

As a second practical application, we will now look at two examples of how 
the curve closing algorithm can be applied to motion data to create periodic 
animations. A similar scheme was developed in [15], but based on an Euler angle 
parametrization for the character joints. For animations exhibiting a large range 
of motions such as rolls and flips, this can result in highly degenerate results 
with characters seemingly just floating in the air. The Lie group based algorithm 
developed in Section 4 exhibits much stronger stability to such outliers. 

In Figure 5, we start with a handspring motion that we repeat three times. 
The plotted curve in the figure shows the trajectory of the right foot of the 
character. Following this trajectory in the upper part of the figure, we can see 
a discontinuity in the foot position when the animation repeats. The problem 
is that the start and end poses of the animation are too different, which results 
in a noticeable jerk when we repeat the handspring. Note that while this gap 
may seem small in the static picture, it is much more noticeable when looking 
at the actual animation. Figure 6 shows the discontinuity in more detail. 

On the bottom of Figure 5, we see the result of applying a few iterations of 
the curve closing method to the handspring animation. The discontinuity has 
been strongly reduced, while the rest of the motion has been preserved. The 
result is an aesthetically much more pleasing motion. 

Figure 7 shows another example of the animation closing method. In this 
case, the character performs a cartwheel, and the animation starts and ends in 
two very different poses, which causes a big discontinuity, particularly in the 
left hand, when the animation is repeated. The closing algorithm manages to 
close this gap to a large extent, while preserving the general appearance of the 
animation. Note that the gap could be closed further by running more iterations 


28 




t 


Figure 5: Application of closing algorithm to a handspring animation. The 
motion is repeated twice. Note how, when repeating the original animation, 
the feet jump backwards to cover a gap between the first and last frame in 
the animation. In the bottom plot, the curve closing method has been used to 
remove this discontinuity. 
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Figure 6: A closer look at the discontinuities in the handspring animation of 
Figure 5. The (blue) left half of both figures shows the last few animation 
sampling points of the animation, the (orange) right half shows the first few 
sampling points when the animation is repeated. To avoid cluttering the figure, 
only the left half of the skeleton has been plotted. 


of the closing algorithm, but at the expense of introducing visual artifacts, such 
as the left hand sliding on the ground during the cartwheel. 

6 Conclusion 

In this article, we have formulated a shape analysis framework for curves on 
Lie groups based on the SRVT approach [42]. This has allowed us to construct 
efficient algorithms to solve two very different problems in computer animation: 
Interpolating between animations and generating cyclic animations. Potential 
further applications include classification and search of animations. 

Future work in this area could involve both joint-wise constraints (e.g., knees 
are not allowed to bend backwards) as well as blending of multiple animations 
using a Karcher mean approach, as described in Section 5.3. Also, investigating 
higher order continuity in the curve closing algorithm could lead to improved 
practical results. 
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Original 
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Figure 7: Application of closing algorithm to a cartwheel animation. Note that 
large different between start and end poses, on the right and the left respectively. 
The motion is repeated once and suffers from a strong jerk when it repeats, 
especially in the left hand. In the second row, the curve closing method has 
been used to alleviate this discontinuity. The remaining gap could be closed 
further using additional iterations of the closing algorithm, but this eventually 
introduces visual artifacts, with the left hand sliding on the ground during the 
cartwheel. 


A Detailed proofs for Section 3 


A.l Proof of Proposition 3.14 

In Proposition 3.14 a formula for the tangent map of the right-logarithnric 
derivative was given. For the readers convenience we repeat the formulation 
of the proposition now. 

Proposition Let c £ C°°(I,G) and v £ T c C°°(I,G) = {w £ C°°(I,TG) \ 
w(t) £ T c i t )G, 'it £ I}. Then writing [•, •] for the Lie bracket in g we have 


Tc rm) = £ t (R£ t) y(v(t))) 


R( C \ t)r (v(t)),s r (cm 


Finally, let us write out the above formula for T c S r in the case that G is a linear 
Lie group: 
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Remark A.l (Proposition 3.14 for linear Lie groups). Let G be a (finite¬ 
dimensional) linear Lie group, c: I —>• G smooth and v £ T c C°°(I,G), i.e. 
v. I —¥ TG with v(t) £ T c ( t )G. Then 


T c 5 r (v)(t ) 


• c(t)) + 


v(t) ■ c 1 (t), c(t) ■ c : (t) 


where products are matrix products and the Lie bracket is the commutator 
bracket. 


Proof of Proposition 3.1 f. Recall from [27, 38.1 Lemma] the product rule for 
the logarithmic derivative: 


$ r (f ' 9)(t) = S r (f)(t) + Ad( f(t))(5 r (g)(t)) for f,g: I G smooth 


Here Ad(y) = (L g )* o (R g - i)* is the adjoint action of G on 0 . Observe that 
the derivative of Ad in the unit e £ G is the derived representation of G on g, 
which is T e Ad(x)(y) = [a;, y], i.e. the Lie bracket in g. Moreover, recall from [17, 
Lemma 2.1] that T eG 5 r (f)(t) = gj£(t), where ec '■ I —> G is the unit in G°°(/, G) 
(cf. Proposition 3.2). Now fix a curve 7 : ] — a, a[xl —► G with 7 ( 0 , t) = c(t ) for 


all t £ / and 
as follows: 


_d_ 

de 


£ — 0 


7 ( 0 1 ) 


v(t). Then the tangent map of S r can be computed 


Tc6 r (v)(t) = - 

d_ 

de 


£—0 


£—0 


FMt,t))= g~ e 


<5 r (7(e,t) • (c(t)) -c(f)) 


£—0 


(t r (y(£,t)-(c(t)) J ) + Ad( 7 (£,t)-(c(t)) *)(d r c 


= s 


£ —0 


(7 (£,*)■ (c(t» 1 


+T - Ad u 


£ —0 


( 7(0 i) • (c(t)) ') ) (<Tc(t)) 


at 


( R (it))*( v (t))) + 


R (c \ t)r (v(t)),6 r (c)(t) 


□ 


( 21 ) 


A.2 Proof of Theorem 3.16 

The aim of this section is to prove Theorem 3.16. Before we do this, we need 
the following auxiliary result. 

Lemma A.2. Let H be a Hilbert space with dim H > 2 and Sh '■= {v £ H \ 
||w|| = 1} its unit sphere. For i = 1,2 fix smooth curves Ui- I —> Sh on a 
compact interval /. 9 Then there is an open subset O C Sh which contains the 
image of u± and U 2 together with a diffeomorphism $: O —> V to a vector space 
V with dim V > 1. 

Proof. We have to distinguish two cases: 

^Recall that the unit sphere of a Hilbert space is a closed submanifold by [32, p. 29 
Example], whence it makes sense to consider smooth curves to Sh- 
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Case 1: 2 <dim H < oo. Since I is an interval and dim Sh = d > 1, 
Sards Theorem (see [32, XVI, §1 Theorem 1.4]) implies that the smooth curves 
ui,U 2 - I —> Sh are not surjective. Hence we can pick x £ Sh such that x £ 
U\(I) U 11 - 2(1 )■ Set O = Sh \ 1®} and let <E>: O —> be the stereographic 
projection through x. Then $ is a diffeomorphism as needed. 

Case 2: H is infinite-dimensional. As H is an infinite-dimensional Hilbert 
space, it is well known (cf. [13]) that Sh is diffeomorphic to H itself. Hence set 
O = Sh and let $: O —> H be the diffeomorphism constructed in ibid. □ 


Notice that the proof of Theorem 3.16 below does not generalise to the case 
dim g = 2. 

Theorem If dim g > 2 then the geodesic distance of C°° (/, g\ {0}) is 
globally given by the L 2 -distance. 


Proof of Theorem 3.16. Let ( 71,(72 £ C°°(/,g\{ 0}). Denote by dc°°(/. £ i\{o}) the 
geodesic distance in C°°(I,g \ {0}). Observe first that the geodesic distance 
coincides locally with the L 2 -distance d^ 2 . In particular, this implies 


dc°°(L g\{o})(9i, 92 ) > d L 2 ( qil q 2 ). 


We have to make sure that both distances coincide even if the minimizing 
geodesic c gii92 (s) = (1 — s)q\ + sq 2 in C°°(I,g) is not contained in C°°(/,g \ 
{0}). Obviously, the minimizing geodesic is not contained in C°°(/,g \ {0}) 
if there are s,f such that c gij92 (s)(f) = 0. Such (s,t) can exist if and only if 


gi (*) 


12W 


i —tttt — —mr* Our aim is now to find a smooth perturbation of ui which 

is arbitrarily close to u\ (with respect to the L 2 -norm) such that the linear 
paths c„ 1>9l (s) = sv 1 + (1 — s)q± and c Vliq2 (s) = sv± + (1 — s)q 2 are contained in 
C°°(/,g \ {0}). We split the problem into two distinct steps: 


Step 1: Construct a perturbation V\ ofqi such that c Vljq2 is in C°° (I , g\{0}) . 
Set u\ := ||~rr and u 2 ■= — to obtain maps which take their image in the unit 
sphere S B = {v £ g | ||u|| = 1}. Since g is a Hilbert space the unit sphere S g is a 
closed submanifold of g (cf. [32, p. 29 Example]) and we see that m £ C°°(I , S g ) 
for i = 1,2. Now it suffices to construct a smooth perturbation v: I —)■ S g of 
u\ such that v(t) 7 ^ u 2 (t) for all t £ I. By Lemma A.2, there is an open set 
O C S g which contains the images of u± and u 2 together with a diffeomorphism 
H>: O —> V to some vector space V with dimV > 1. Consider now the smooth 
curve w\ = $ o u\ — •!> o u 2 '■ I —> V. Since dim V > 1 we can clearly construct 
a smooth mapping w: I —> V \ {0} which satisfies 


sup||$ o u\ — (w + $ o u 2 )|| = sup||u;i(t) — w(t)\\ < 6 (22) 

tei tei 


for some arbitrary but fixed S > 0 (the control S will be needed in Step 2 below). 
Then v := <k -1 o (w + $ o u 2 ) : I —> S g is smooth, and satisfies v(t) ^ u 2 (t) for 
all t £ I (since <k(u) = w(t) + 4> ou 2 (t) 7 ^ >I > o'U 2 (t))- Define V\(t) := v(t)-\\qiW ( t ) 
to obtain a smooth map such that c Vl q2 (s ) = sv 1 + (1 — s)u 2 is a path in 

C°°(I,Q \ {0}). 

Step 2: Adjust v\ such that also c Vim is in C°°(I , g \ {0}). Since q ±: I —> 
g \ {0} is smooth and I is compact, r := inf te /||( 7 1 (t)|| > 0. Hence if 

sup||( 7 i(f) - vi(t)\\ < r, (23) 

tei 
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then c„ li9l is contained in C°°(I,g \ {0}). To see that we can choose v\ = 
||gi|| • + $ o u 2 ) such that (23) is satisfied, we consider the map 

e-. C(I, V ) -+ C(I ,0 \ { 0 }), / ^ (t ^ ||gi(i)|| (*-\m + $ o U 2 (t)))) 

and from Step 1 we have v\ = 9(w). We claim that 9 is continuous if we endow 
the spaces of continuous maps C(I,V) and C(I,g \ {0}) with the topology 
induced by WfW^ '■= sup t 6 / ||/(f)|| (where ||-|| denotes the norm of V and g, 
respectively). If this is true, then the proof can be completed as follows. Since 8 
is continuous with = 9 i> we can choose S r > 0 so small to ensure 

that (23) holds if (22) is satisfied. Thus c vi<qi takes its image in C°°(I,g\ {0}) 
if S r > 0 is small enough. 

Summing up, we have seen that we can always construct a smooth pertur¬ 
bation vi of qi such that the linear paths c Vl , qi and c v , m are contained in 
C°°(I, g \ {0}). Moreover, for e > 0 the estimates in Step 2 show that we can 
choose 5 r > S e > 0 such that the left hand side of (23) is smaller than e. Hence 
for each e > 0 we can choose a smooth map i>f such that di 2 (uf,gi) < e. In 
particular, v\ converges to q\ with respect to the L 2 -distance. Then the geodesic 
distance satisfies 

dc“(/, 0 \{o})( 9 i> 82) < d c *>(i, g \{o})(qi,vl) + d c °°(i, g \{o}){vl,q 2 ) 

= d L 2(q 1 ,vl) + d L 2 (t;f, q 2 ) 

< £ + d L 2(v{,q 2 ) d L 2(q 1 ,q 2 ). 

Thus the geodesic distance of 0 \ {0}) coincides with the L 2 -distance. 

Proof of the claim: 6 is continuous. Recall that the topology on C(I, V ) 
induced by 11 • 11 ^ coincides with the compact open topology. Hence [14, Theorem 
3.4.2] shows that the map 7 $: C(I, V ) —> C(I, S B ),f 4 > -1 o / is continuous 
as is continuous. Further, C°°(I,V) with the above topology is a Banach 
space, whence 


h- C(I,V)C(I, S g ), f e->• 7 $ ° (/ + "h ° iti) 

is continuous. Now as ||||gi||< H^iHooll/Hoo we deduce that n qi : C(I,S g ) —¥ 
0(1, g \ {0}), / i-)-1|gi || / is continuous. In conclusion, 9 = n qi o h is continuous. 

□ 
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